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We use simple analytic considerations and a Monte Carlo calculation of nucleons in a box to 
argue that the use of Jastrow functions as short-range correlators in the commonly employed two- 
body-cluster approximation causes significant errors in the matrix elements for double-beta decay. 
The Jastrow approach appears to agree with others, however, if many-body clusters are included. 
A careful treatment of the charge-changing analog of the nuclear pair density shows, in addition, 
that differences between Unitary Correlator Operator Method and Brueckner methods for treating 
short-range correlations in double-beta are less significant than suggested by previous work. 



I. INTRODUCTION 

New experiments to measure the rate of neutrinoless 
double-beta {01^/3(3) decay will provide information about 
neutrino masses if neutrinos are Majorana particles. Un- 
fortunately, one must know the value of the nuclear ma- 
trix element governing the decay to extract that informa- 
tion from an observed rate (or rate limit) For that 
reason, attempts to better calculate the matrix elements 
appear regularly in the literature. 

The matrix elements involve products of one-body de- 
cay operators and a sum over intermediate states, but 
the closure approximation allows them to be represented 
to high accuracy 0] by the ground-statc-to-ground-state 
transition matrix element of a two-body operator 



(1) 



a<b 



The matrix elements Mfi = (/| Ai'^'^ \i) can therefore be 
affected by the strong repulsive correlations that alter 
pair distribution functions at short distances. 

For many years, theorists were satisfied to simulate 
these correlations by using a phenomenological Jastrow 
function f(rab) in the two-body cluster approximation to 
modify the operator A^: 



Mab =^ f{rab)Mabf{rab) 



(2) 



where Tab = |ra — r;,| is the magnitude of the distance 
between the two nucleons. The Jastrow function almost 
always had the form prescribed in Ref. [1] : 



(3) 



with r in fm. Recent work has questioned this 
prescription. Refs. treated short-range correla- 

tions through the Unitary-Correlation Operator Method 
(UCOM), which has the advantages of wave- function 



overlap preservation and a range of successful applica- 
tions Refs. [1| and Q computed the effects of short- 
range correlations within well-defined Brueckner-based 
approximation schemes. All these papers found smaller 
effects on matrix elements than the phenomenological 
Jastrow function yields. Because they all limited their 
analysis to two-body correlations, however, their predic- 
tions for the size of short-range effects do not come with 
an iron clad guarantee. 

In fact, all these methods imply the existence of many- 
body effects that are always neglected in applications 
to double-beta decay. In Jastrow-based treatments, our 
subject here, the full correlated wave function is written 
schematically as 



n ^ 

\a<b<c 
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where is a Slater determinant or a generalization 
thereof, Fab is a two-body Jastrow correlator depend- 
ing on rab and the spins and isospins of particles a and 
6, and Tabc is a similar three-body correlator, which we 
will ignore from here on. In recent years, practitioners 
have developed a range of techniques for moving beyond 
the two-body cluster approximation in Eq. ^ and in- 
cluding many-body correlations generated by the prod- 
uct of F's (in addition to explicit three-body correlations 
generated by a single T) in Eq. Q. Cluster expansions 
and the Fermi-hypernetted-chain method (see, e.g., Refs. 
[lol ITl| and references therein) include three-and-more- 
body clusters, and quantum Monte-Carlo methods allow 
an evaluation of the contributions of all clusters. The 
Jastrow approaches now yield accurate observables, in- 
cluding two-body density distributions with short-range 
correlations, in light nuclei [l^ and nuclear matter [l3j . 
Here, after analyzing the two-body cluster approxima- 
tion in /3/3-decay, we see whether an initial application of 
quantum Monte Carlo with many-body correlations in- 
cluded supports the phenomenological two-body-cluster 
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Jastrow method used traditionally, or whether it sup- 
ports one or more of the approaches introduced recently. 
We also point out that apparent differences between the 
results of UCOM and Brueckner methods are largely fic- 
titious. 

Heavy nuclei are still too complicated for Monte-Carlo 
methods in their current forms, so to evaluate many-body 
Jastrow effects we look instead at a simplified version of 
asymmetric nuclear matter. We make this choice with 
the idea that short-range correlations are nearly univer- 
sal in nature, depending little on longer-range structure 
of the environment in which the correlated nucleons are 
embedded, provided that environment has the correct 
density. 



II. TWO-BODY CLUSTER APPROXIMATION 

In the S ~ Q T ~ 1 channel that determines the con- 
tribution of short distances to the /3/3 amplitude, realis- 
tic variationally-detcrmined correlation functions Fab arc 
not so different from the Miller-Spencer Jastrow function. 
Figure [1] shows a typical nuclear- matter example, follow- 
ing the calculations of Ref. [l^, alongside the Miller- 
Spencer function and the effective scaling function, ob- 
tained from the ratio of calculations with and without 
short-range correlations, that appears in the Brueckner- 
based work of Ref. Q All the functions go to unity at 
large r, but the Brueckner-based function has a sizeable 
"overshoot" near ?■ = I fm. The Miller-Spencer func- 
tion has a much smaller overshoot (occurring at larger r, 
which is made less important by the radial falloff of the 
Oi'pf] operator) leading to a significantly smaller 01/(5(3- 
matrix element. The variational nuclear- matter resem- 
bles the Miller-Spencer function but has essentially no 
overshoot, and so if applied like that function via Eq. ^ 
it will produce an even smaller matrix element. 

The use of the F from Eq. (|4]) to multiply a two-body 
operator as in Eq. ^ is often called the two-body clus- 
ter approximation, because all terms arc discarded except 
those in which the transition operator and the correlators 
act on the same pair of particles. This approximation ap- 
pears to be reasonably good for number-conserving two- 
body densities. The dot-dashed line in Fig. [T] displays the 
distribution goi(^) in the 5 = 0,T = 1 channel, follow- 
ing Ref. [ll] , which incorporates the ful product over all 
pair correlations of Eq. (jH) . This full goi (r) is somewhat 
smaller than the corresponding F'^ because many-body 
tensor correlations promote a fraction of the spin-singlet 
pairs to spin-triplet pairs, so that the number of singlet 
pairs is reduced. The reduction has also been seen in 
light nuclei [13] , though the corrections are not large ei- 
ther there or here. 

In (3(3 decay, the picture must be different, however. 
To see why, consider the charge-changing analog of the 



(spin-independent) two-body density: 

P^(r)^(/|^<5(r-r,,)r+r+|») , (5) 

a<b 

where F stands for Fermi. If we weight this function 
with Hpir), the radial part of the Fermi 0v(3(3 operator 
(given approximately by l/r), and integrate, we get the 
Fermi piece of the 0v(3(3 matrix element. If we integrate 
Ppir) without any weighting, we get (/| Y^aKb^tTb N)' 
which must vanish because the isospins of \i) and |/) are 
different (in the very good approximation that isospin is 
conserved), while the operator between them is propor- 
tional to the square of the isospin-raising operator. 

Figure [2] shows Pp (r) for the shell- model calculation 
of the /3/3-decay of ^^Se in Refs. [H and [Hi. The sohd 
curve contains no Jastrow function and has area of zero 
beneath it. The dashed curve is the result of of the 
Brueckner-based calculations in Ref. Its overshoot 
at r just greater than one causes the integral to stay 
very close to zero despite the suppression at very small 
r. But the use of the two-body Jastrow function Fqi 
a la Ref. [l^l (dotted curve) suppresses contributions at 
small r without an overshoot and thus leads to an in- 
tegral of 0.006. Substituting the pair distribution func- 
tion goi would only make the problem here worse. The 
Millcr-Spcnccr Jastrow function yields a little bit of over- 
shoot but not nearly enough, and results in an integral 
of 0.0075. 

It seems, then, that a realistic treatment of short-range 
correlations must yield an overshoot in Ppif) if it is to 
preserve isospin (The UCOM procedure does this exactly, 
by construction). When Jastrow functions arc extended 
beyond the two-body cluster approximation, the effec- 
tive functions that result must therefore look different 
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FIG. 1. (Color online) Squares of Jastrow functions Fat from 
calculations following Ref. [l^l (dotted black line, spin-singlet 
only), from Miller and Spencer [3l (solid red line) and from a 
fit to the results of a microscopic Brueckner-based calculation 
[^ (dashed blue line) . The purple dot-dashed line comes from 
three- and more-body corrections to the dotted line. 
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FIG. 2. (Color online) The charge-changing two-body de nsity 
Ppir) for the shell-model calculation of ®^Se in Refs. p^. [lq| . 
The solid red line is the result without short-range correla- 
tions, the dashed blue line is that from the Brueckner-based 
calculation of Ref. @] the dotted black line applies the Jas- 
trow function from the approach in Ref. [T^ in the two-body- 
cluster approximation, and the dot-dashed purple line applies 
the Miller-Spencer Jastrow function. The inset magnifies the 
left-hand part of the figure. 



for charge-changing densities, which involve only valence 
nucleons, than for like-particle densities, to which all nu- 
cleons contribute coherently. 

It is not hard to get an idea of how this happens. Let us 
consider spin-and-isospin-independent two-body correla- 
tors Fab (with no three-body correlators Tabc) in Eq. (|4|) 
and a general charge-changing operator Mat- Writing 
F^j^ = 1 -|- hab, and expanding to first order in h gives 

{f\M\t) = {.fo\Y,Mabl[{l + h,d)\^o) (6) 

a<ib c<d 

ab 1^0/ 

a<6 a<ib 

+ (./o| ^ MabiKc + hbc) 

a<b 

c^{a,fc} 

-f(/o| Mabhcd\io) + O{h^) 

a<b 

c<d^a,b 

-(/ol^A^ab U+Y.^-d) h)+Oih^) 
a<b \ c<d / 

where |io) and |/o) are Slater determinants, and in the 
third and fourth lines, hmn = h„rn if n < m. 

The second line in Eq. ([6]) involves the bare two-body 
transition operator and the two-body-cluster correction. 
The third line involves an effective three-body opera- 
tor, and the fourth line an effective four-body operator. 
Terms of higher-order in h generate even higher many- 
body operators. 



Now let the neutron number exceed the proton number 
so that the Slater determinants |io) and |/o) have well- 
defined isospins that differ from each other, and consider 
the operator A4ab = '''^'''b' ■ The matrix element above is 
then the integral of Ppir), i.e. zero. Although the first 
term in the second line indeed gives zero, the second, as 
we've seen, does not. The inclusion of all terms first- 
order in h must restore the value zero, however, because, 
as the last line shows, the result can be obtained by act- 
ing on |zo) with the isospin- preserving two-body operator 
Sa<h ^"fc before acting with the transition operator. It 
is not hard to show that effective four-body term con- 
tributes the same amount as the two-body-cluster correc- 
tion, and the effective three-body term contributes twice 
that amount with the opposite sign, so that the sum of 
terms indeed vanishes. But this also means that, at least 
to first order in h, three- and-four body effective oper- 
ators are just as important for the quantity J PF{r)dr 
as is the effective two-body correction generated by the 
two-body cluster approximation. This perhaps surpris- 
ing conclusion leads us to examine the charge- changing 
density itself and the higher order corrections in a model 
amenable to numerical solution. 



III. MANY NUCLEONS IN A BOX 

We consider a cubic box with each side L = 4.85 fm 
and periodic boundary conditions. In the box are 2 pro- 
tons and 16 neutrons (so that the nucleon density is very 
near nuclear- matter density), which decay to 4 protons 
and 14 neutrons. The protons in the initial state and 
all but the last two neutrons in that state are in filled 
fermi levels, and the last two (valence) neutrons are in 
the spin-zero two-body pairing wave function: 

|^„)=AA Y |k,-k;^ = 0), (7) 

where v stands for valence. A/" is a normalization con- 
stant, and the set K contains vectors in which two k 
components arc equal to ±2t: / L and the third is zero. 
In the final state the neutrons and all but the last two 
protons are in filled fermi levels; the two valence protons 
are in the configuration above, but with the set K 
containing vectors with one component equal to ±211 /L 
and the other two equal to zero. 

We use quantum Monte Carlo to evaluate Ppir) be- 
tween states of the form Eq. (|4|), where now the states 
|zo) and |/o) are those just described, and with simple 
spin-and-isospin two-body correlators Fab and no three- 
body correlators Tabc, as considered previously. Figure [3] 
shows the results with the Miller-Spencer Jastrow func- 
tion. The two-body-cluster approximation again has very 
little overshoot, but the full results, including clusters of 
all size, has considerably more, so that the integral van- 
ishes as it should. Also shown is the result with the 
effective Jastrow function fit to the Brueckner-based cal- 
culation of Ref. Q (which was done in finite nuclei). It 
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FIG. 3. (Color online) Monte-Carlo calculation of Ppir) for 
2 protons and 16 neutrons in a box decaying to 4 protons and 
14 neutrons. The red circles are the result with no Jastrow 
correlators, the purple diamonds include the Miller-Spencer 
Jastrow correlator in the two-body-cluster approximation, the 
black upward-pointing triangles are the full result with that 
correlator, and the blue downward-pointing triangles apply 
the effective Brueckner based two-body-cluster Jastrow from 
Ref. m. The numerical error associated with the values of 
Ppir) are usually smaller than the size of the corresponding 
symbols. 



is now quite close to the full many-body Miller-Spencer 
result, and the remaining discrepancy is probably mostly 
due to the simplicity of our model. Surprisingly, the use 
of a Jastrow function with no overshoot at all gives al- 
most the same result as the Miller-Spencer function when 
many-body correlations are taken fully into account. 

To show the effects of these various functions on 
decay, we define functions Ck[t), K = F,GT, for the 
Fermi and Gamow- Teller parts of the Oiz/Sf] operator. (If 
we wanted to be accurate we would also define one for the 
very small tensor term.) These functions are the products 
of the densities Ppif) and the analogous density 

PgtW ^ (/I ^ <5 (r - r,b) a, • a,r+r+ \i) (8) 

a<b 

with functions Hp{r) and Hgt{t) that specify the radial 
dependence of the operators. In other words 



CK{r) ^ HK{r)PK{r) ., 



(9) 



with 



HF{r) = HoTir) 



2R 



dq 



sm qr 



q + E-{E,+Ef)/2 

(10) 



The quantities E, Ei, and Ef are energies to which the 
Hk are not very sensitive. Equation (fTUI) holds only if we 
neglect nucleon form factors and forbidden currents; in 
the complete, more complicated expressions Hp ^ Hqt 



|17| . We use the simplified forms here because they are 
sufficient to make our point. 

Figure [4] displays Cpir) from the Monte Carlo. (In 
this simple calculation Cqt is just proportional to Cp 
because the correlator is spin-independent and the two 
valence nucleons that participate in the decay are locked 
into a spin-zero configuration). The full solution clearly 
corrects the extreme suppression of the matrix ele- 

ment created by the two-body-cluster approximation, in 
a way consistent with the analysis of the integral in sec- 
tion ini Differences with the Brueckner treatment are 
fairly small and due once again at least in part to the 
unusual system we analyze here. Effects beyond the two- 
body cluster approximation are thus both required and 
apparently sufficient to describe short-range correlations 
in /3/3 decay. 




FIG. 4. (Color online) The curves Crir) corresponding to 
the distributions Ppir) of Fig. [S] without forbidden currents 
or nucleon form factors. The symbols indicate the same ap- 
proximations as in Fig. 3. 



IV. UCOM 

In this section we digress from our main line of inquiry 
to take up apparent differences between the results of 
Brueckner methods and UCOM. Fig. 2 of [3], Fig. 9 of 
Ref. [Ol and Fig. 4 of 0] present Oi^/3/3 distribution func- 
tions (similar to the function ir) presented here in 
Fig. [4]) with the UCOM (and other) treatments of short- 
range correlations. Unlike the Brueckner-based curves in 
our Figs, m ~ [3] the UCOM curves show no overshoot. 
But the reason is that the "contribution from distance 
r" has been treated differently when UCOM correlations 
are considered than when other methods are used. The 
UCOM prescription requires that the operator Tab be re- 
placed by a shifted version R+{rab) (where the function 
-R-l- is usually determined variationally) in any opera- 
tor that doesn't depend on momentum. Thus, to eval- 
uate the Qv^P matrix element, one replaces Hxirab) by 
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HK{R+{rab)). 

The shifting imphes that UCOM produces a function 
Cf (which we use as an example because it is simpler 
than Cgt) given by 

Cj^ir) = E (/I HF{R+irabmr - R+ir,b))T+T+ \z) 

a<b 

= Hpir) (/I J2 - R+irab))T+T+ \z) 

a<b 

^HAr)P^{r), (11) 

(where U stands for UCOM). In prior work on UCOM 
/3/3-decay, however, the "(Fermi) contribution from a 
given r" was defined instead by simply replacing Tab with 
R+{rab) in H, viz: 

C%'{r) = (/I E Hp{R+{r,,bmr - rab)r+T+ |z) 

a<b 

^ HF{R+{r))PF{r) . (12) 

This definition, which leaves rab unshifted in the delta 
function, gives the correct result for the Fermi matrix 
element when r is integrated over, but does not define 
an observable and is not what is calculated in other ap- 
proaches. The correct expression, Eq. (fTTj) . is a bit more 
complicated to evaluate, but has a pronounced overshoot. 
Figure[5]comparcs the distribution CoTir) from the shell- 
model-Brueckner treatment of ®^Se — > ^^Kr in Ref. Q 
(with no forbidden currents or nucleon form factors) to 
the properly defined UCOM distribution for the same de- 
cay. The two curves are extremely close to one another, 
and quite different from CQ^{r), which is also shown. 
The UCOM and Brueckner pictures are therefore more 
similar than previously thought. 
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FIG. 5. (Color online) The function CcTir), defined in Eqs. 
(IH)) and solid red curve), the corresponding UCOM func- 
tion CoT{r), defined as in Eq. (dashed blue curve), the 
previously-used UCOM function Cg^(r), defined as in Eq. 
(|12|l (dot-dashed purple curve), and the Brueckner-based ver- 
sion (dotted black curve). The new UCOM and Brueckner 
curves are very similar. 



V. DISCUSSION 



The main point of this paper, to which we now re- 
turn, is that the use of Jastrow functions in the two- 
body-cluster approximation suppresses short-range con- 
tributions too much, and that the problem is fixed by 
including many-body correlations. This discovery raises 
the question of whether existing treatments are adequate. 
They include long-range many-body correlations in shell- 
model or QRPA wave functions but allow only two parti- 
cles to be correlated at short distances. Is that sufficient? 

It is hard to answer the question definitively because 
the approach taken here is so different from the others. 
We can say that the very-short-range correlations arc un- 
likely to be altered; our figs. [3] and [4] show that correc- 
tions to the two-body-cluster approximation are barely 
noticeable below about r = 0.7 fm. But corrections are 
large at 1 fm or so. It is far from obvious that the 
marriage of UCOM or Brueckner treatments of short- 
range correlations to shell-model or QRPA treatments of 
longer-range correlations incorporates all important ef- 
fects at the intermediate range r « 1 fm. The UCOM 
procedure generates three- and higher-body correlations 
that have been neglected in almost all applications to 
date, and the Brueckner-based double-beta work has so 
far not included contributions from, e.g., three-particle 
ladders. As for the shell model and QRPA, they leave 
untreated a significant range of single-particle energies 
between those contained in the calculation and those rep- 
resented as short-range effects. Whether these omissions 
are significant is still an open question. 

In the meantime, however, it appears that the UCOM 
and Brueckner methods give reasonable short-range cor- 
relations. Like the full Jastrow calculations here, they 
supply the overshoot required to preserve isospin sym- 
metry. Higher-body corrections in these schemes appear 
unlikely to be as large as they are in the Jastrow ap- 
proach, which violates isospin symmetry in the two-body 
cluster approximation. Short-range effects in /3/3-decay 
thus seem to be mostly under control. 
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